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Abstract 

We consider multi-loop integrals in dimensional regularisation and the corresponding Lau- 
rent series. We study the integral in the Euclidean region and where all ratios of invariants 
and masses have rational values. We prove that in this case all coefficients of the Laurent 
series are periods. 



1 Introduction 



The calculation of loop integrals in perturbative quantum field theories is essential for accu- 
rate and precise theoretical predictions to be compared to experiments. The calculation of loop 
integrals is complicated by the occurrence of ultraviolet and infrared singularities. Ultraviolet 
divergences are related to the high-energy behaviour of the integrals. Renormalisation absorbs 
these divergences into renormalisation constants. Infrared divergences may occur if massless 
particles are present in the theory. For infrared-safe observables they cancel in the final result, 
when summed over all degenerate states. 

Dimensional regularisation [1-3] is usually employed to regularise these singularities. Within 
dimensional regularisation one considers the loop integral in D space-time dimensions instead 
of the usual four space-time dimensions. The result is expanded as a Laurent series in the 
parameter e = (4 — D)/2, describing the deviation of the D-dimensional space from the usual 
four-dimensional space. The singularities manifest themselves as poles in 1/e. Each loop can 
contribute a factor 1/e from the ultraviolet divergence and a factor 1/e 2 from the infrared di- 
vergences. Therefore an integral corresponding to a graph with / loops can have poles up to 
l/e 2/ . 

In this paper we consider multi-loop integrals in dimensional regularisation and the corre- 
sponding Laurent series. It is an interesting question to ask, what type of numbers the coeffi- 
cients of the Laurent series are. These numbers certainly depend on the momenta of the external 
particles and the masses of the particles propagating in the loops. In order to keep the loop inte- 
gral dimensionless, one introduces an additional arbitrary mass /u and the final result will depend 
on this mass as well. For scalar loop integrals the dependence on the external momenta is only 
through Lorentz invariants like (pi + P2 + ■■■ + Pj) 2 - Let us assume that all ratios of these in- 
variants and masses are rational numbers. They are certainly real numbers. Since the rational 
numbers are dense in the real numbers, this is no severe restriction. Let us further assume that 
we consider the loop integral in the Euclidean region, meaning that all masses are positive or 
zero and all invariants are negative or zero. From explicit calculations we know, that for all one- 
loop integrals the e°-coefficient involves only rational numbers, logarithms and dilogarithms. In 
two-loop integrals we encounter multiple poly logarithms, which have been studied extensively in 
the literature by physicists [4-17] and mathematicians [18-32]. For the massless two-loop two- 
point function it has been shown that all coefficients of the Laurent expansion can be expressed in 
terms of multiple zeta values [33]. As a last example let us mention that there is strong numerical 
evidence that the three-loop two-point function with equal internal masses involves elliptic inte- 
grals [34]. We are therefore tempted to ask, what is the unifying theme for all these examples? 
Is there a common set of numbers or functions, to which all these examples belong? We observe 
that common to all examples is the fact, that these functions yield periods when evaluated with 
rational numbers as arguments. Periods are special numbers and are introduced in section[2] This 
raises the question if the observation we made for the few examples above can be shown to hold 
for all loop integrals. 

In this paper we show that this is indeed the case and prove a rather general result: With the 
assumptions above (Euclidean region and all ratios of invariants and masses rational) we show 
that all coefficients of the Laurent series of an arbitrary scalar multi-loop integral are periods. 
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This result puts strong restrictions on the class of functions, which can appear in the calcula- 
tion of multi-loop integrals. It is well known that a general multi-loop integral involving tensor 
structures in the numerator can be reduced to scalar integrals. An algorithm for this reduction is 
reviewed in section|4j Therefore our result extends to all multi-loop integrals. 

Our result extends and generalises a theorem by Belkale and Brosnan [35] on the Laurent 
expansion of Igusa local zeta functions. Our proof is constructive and based on recent work on 
the resolution of singularities for multi-loop integrals [36]. The proof uses the close relationship 
of Feynman integrals with algebraic geometry, in particular the equivalence of iterated sector 
decomposition [37-40] withHironaka's polyhedra game [41,42]. Methods of algebraic geometry 
have been introduced in the context of Feynman integrals in [43]. 

This paper is organised as follows: In section [2] we introduce the set of periods. Section [3] 
reviews basic facts about scalar multi-loop integrals. Section @] shows that all loop integrals 
occurring in a quantum field theory can be reduced to these scalar integrals. In section [5] we state 
the main theorem of this paper and present a proof. Section [6] contains the conclusions. In the 
appendices we discuss Hironaka's polyhedra game and the Laurent expansion of the prefactors 
accompanying the loop integrals. A third appendix discusses the relation to the work of [43] for 
the sequence of blow-ups. 



Periods are special numbers. Before we give the definition, let us start with some sets of numbers: 
The natural numbers N, the integer numbers Z, the rational numbers Q, the real numbers M 
and the complex numbers C are all well-known. More refined is already the set of algebraic 
numbers, denoted by Q. An algebraic number is a solution of a polynomial equation with rational 
coefficients: 



As all such solutions lie in C, the set of algebraic numbers Q is a sub-set of the complex numbers 
C. Numbers which are not algebraic are called transcendental. The sets N, Z, Q and Q are 
countable, whereas the sets R, C and the set of transcendental numbers are uncountable. 

Periods are a countable set of numbers, lying between Q and C. There are several equivalent 
definitions for periods. Kontsevich and Zagier gave the following definition [44]: A period is a 
complex number whose real and imaginary parts are values of absolutely convergent integrals of 
rational functions with rational coefficients, over domains in M." given by polynomial inequalities 
with rational coefficients. Domains defined by polynomial inequalities with rational coefficients 
are called semi- algebraic sets. 

We denote the set of periods by P. The algebraic numbers are contained in the set of periods: 
Q G P. In addition, P contains transcendental numbers, an example for such a number is 7t: 



2 Periods 
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The integral on the r.h.s. clearly shows that n is a period. On the other hand, it is conjectured 
that the basis of the natural logarithm e and Euler's constant Je are not periods. 

We need a few basic properties of periods: The set of periods P is a Q-algebra [44,45]. In 
particular the sum and the product of two periods are again periods. This is immediate for the 
multiplication of two periods: Let a and b be two periods, given through 

a = J d n xf(x), b = f d m yg(y), Gi C R M , G 2 C R m , (3) 

G\ G2 

where f(x) and g(y) are rational functions with rational coefficients. Then 

a-b = J d n xd m y[f(x)g(y)]. (4) 

G x xG 2 

For the sum of two periods we consider 

G = Gi x{0}x [0,l] m U [0,l]"x{l}xC 2 C l"xlxr. (5) 

Then 

a + b = Jd n xd m y[(\-t)f(x)+tg(y)], (6) 

G 

where t is the coordinate of the middle factor lofR"xRx R m . 

The defining integrals of periods have integrands, which are rational functions with rational 
coefficients. For our purposes this is too restrictive, as we will encounter logarithms as integrands 
as well. The following lemma is easy to prove: 

Lemma 1: Let G C R n be a semi-algebraic set and f(x) and g(x) two rational functions with 
rational coefficients. Assume that the integral 

/ = Jd n xf(x)lng(x) (7) 

G 

is absolutely convergent. Then / is a period. 
Proof: We have 

g g a 

where G' C R" +1 and (jq, ...,x n ,t) G G' if (xi,...,x n ) e G and t G [0, 1]. Clearly G' is again a 
semi-algebraic set. Therefore the integral / is a period. 

With the same technique of introducing additional variables one shows for rational functions 
fi (x), g\ (x), fi{x) and g2{x), all of them with rational coefficients, the following: If the integral 

J = Jd n x{f 1 (x)lng 1 (x)+f 2 (x)\ng 2 (x)} (9) 

G 
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is absolutely convergent, then it is a period. Here we have 



a 



d"x dti dtl | ^'M-'j + I , (10) 



with (jci,...,jc„,fi,f2) eG'cR" +2 if (xi,...,x n ) GG and t\ G [0,1] andf 2 G [0,1]. This shows that 
G' is a semi-algebraic set and therefore 7 is a period. 

As a final example we consider with the definitions as above the integral 



K = fd n xf{x) lngi(x) \ng 2 (x) 



G 

1 1 



y * *y * y /w i00 _ 1Vi+1 (>a(jc) _ 1)fa+1 - (") 







If the integral K is absolutely convergent, then it is a period. 

Clearly these examples can be combined and iterated and we conclude that in the defining 
integral for a period also integrands, which are linear combinations of products of rational func- 
tions with logarithms of rational functions, all of them with rational coefficients, are allowed. 



3 Feynman integrals 

In this section we introduce the central object of our investigations, Feynman loop integrals. 
Most of the material in this section is well-known to physicists and can be found in many text- 
books [46-49]. This section is included for mathematicians with little background on Feynman 
diagrams, but interested in the most important facts. In this section we restrict ourselves to 
scalar integrals. These are integrals, where the numerator of the integrand is independent of the 
loop momentum. A priori more complicated cases, where the loop momentum appears in the 
numerator might occur. However, there is a general reduction algorithm, which reduces these 
tensor integrals to scalar integrals [50,51]. This reduction algorithm is reviewed in section HI 

To set the scene let us consider a scalar Feynman graph G. Fig. Q] shows an example. In 
this example there are three external lines and six internal lines. The momenta flowing in or out 
through the external lines are labelled pi, pi and pi and can be taken as fixed vectors. They are 
constrained by momentum conservation: If all momenta are taken to flow outwards, momentum 
conservation requires that 

PI+P2+P3 = 0. (12) 

At each vertex of a graph we have again momentum conservation: The sum of all momenta 
flowing into the vertex equals the sum of all momenta flowing out of the vertex. A graph, where 
the external momenta determine uniquely all internal momenta is called a tree graph. It can be 
shown that such a graph does not contain any closed circuit. 
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In contrast, graphs which do contain one or more closed circuits are called loop graphs. If we 
have to specify besides the external momenta in addition / internal momenta in order to determine 
uniquely all internal momenta we say that the graph contains / loops. In this sense, a tree graph 
is a graph with zero loops and the graph in fig. Q] contains two loops. Let us agree that we label 
the / additional internal momenta by k\ to 

Feynman rules allow us to translate a Feynman graph into a mathematical formula. For a 
scalar graph we have to substitute for each internal line j a propagator 

(13) 



2 2 ' 



Here, qj is the momentum flowing through line j. It is a linear combination of the external 
momenta p and the loop momenta k: 

Qj = 1j(p,k). (14) 

nij is the mass of the particle of line j. The Feynman rules tell us also to integrate for each loop 
over the loop momentum: 

d A k 

(15) 



However, there is a complication: If we proceed naively and write down for each loop an inte- 
gral over four-dimensional Minkowski space, we end up with ill-defined integrals, since these 
integrals may contain ultraviolet or infrared divergences. Therefore the first step is to make these 
integrals well-defined by introducing a regulator. There are several possibilities how this can be 
done, but the method of dimensional regularisation [1-3] has almost become a standard, as the 
calculations in this regularisation scheme turn out to be the simplest. Within dimensional regular- 
isation one replaces the four-dimensional integral over the loop momentum by an D-dimensional 
integral, where D is now an additional parameter, which can be a non-integer or even a complex 
number. We consider the result of the integration as a function of D and we are interested in the 
behaviour of this function as D approaches 4. It is common practice to parameterise the deviation 
of D from 4 by 

D = 4-28. (16) 

The divergences in loop integrals will manifest themselves in poles in 1/e. In an /-loop integral 
ultraviolet divergences will lead to poles 1 /e at the worst, whereas infrared divergences can lead 
to poles up to 1/e . We will also encounter integrals, where the dimension is shifted by units of 
two. In these cases we often write 

D = 2m -2e, (17) 

where m is an integer, and we are again interested in the Laurent series in e. 

Let us now consider a generic scalar /-loop integral Iq in D = 2m — 2e dimensions with 
n propagators, corresponding to a graph G. Let us further make a slight generalisation: For 
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Figure 1: An example of a two-loop Feynman graph with three external legs. 



each internal line j the corresponding propagator in the integrand can be raised to a power V/. 
Therefore the integral will depend also on the numbers Vi,...,v„. In this paper it is sufficient to 
consider only the case, where all exponents are natural numbers: V/ G N. We define the Feynman 
integral by 

n 

nr(v,-) lo/2 r'dPkr yi 1 nsl 

with v = Vi + ... + v„. The momenta qj of the propagators are linear combinations of the external 
momenta and the loop momenta. In eq. (fT"8l) there are some overall factors, which are inserted 
for convenience: There is a prefactor consisting of Gamma-functions. The arbitrary mass scale 
H ensures that eq. (fT"8l) is dimensionless. The integral measure is now d D k/{m D l 2 ) instead of 
d D k/ '(2n) D ', and each propagator is multiplied by i. The prefactors are chosen such that after 
Feynman parameterisation the Feynman integral has a simple form. The change of prefactors 
and the Laurent expansion of the neglected prefactors is discussed in detail in the appendix [Bj 

How to perform the D-dimensional loop integrals ? The first step is to convert the products 
of propagators into a sum. This can be done with the Feynman parameter technique. In its full 
generality it is also applicable to cases, where each factor in the denominator is raised to some 
power Vj. The formula reads: 



/ |7 / n r(v y ) 1 



xj>0 




(19) 



Applied to eq. (PTEl) we have 

j^XiPi = j^Xii-qf + mf). (20) 
i=\ i=\ 

One can now use translational invariance of the D-dimensional loop integrals and shift each loop 
momentum k r to complete the square, such that the integrand depends only on k 2 . Then all D- 
dimensional loop integrals can be performed. The master formula for a single loop integration 
is: 



/ 



d D k (-k 2 r F(f +q)r(v-|-q) t/-§-« 

iK%[-uk*+F} v " r(f) r(v) F v-q-a 
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The functions U and F depend on the Feynman parameters and are obtained after Feynman 
parametrisation from completing the square. In eq. (I2TT) we allowed additional powers (— k 2 ) a of 
the loop momentum in the numerator. This is a slight generalisation and will be useful later. 

With eq. (|2Tj) one can perform all D-dimensional loop integrals iteratively. As the integrals 
over the Feynman parameters still remain, this allows us to treat the D-dimensional loop integrals 
for Feynman parameter integrals. One arrives at the following Feynman parameter integral, 
which is the central object of investigation of this article: 

Ig = (M) I dx8(l-^Xi)\\\x/ 1 jV _ /Z)/2 • (22) 



Xj >o 



The functions 11 and J depend on the Feynman parameters. If one expresses 

n III 

Y^Xji-qj + mj) = -££^M r A+£2^-e r -7, (23) 

j= 1 r= 1 s= 1 r= 1 

where M is a / x / matrix with scalar entries and Q is a /-vector with four-vectors as entries, one 
obtains 

U = det(M), 7 = det(M) (-J + QM~ l Q) . (24) 

Alternatively, the functions U and f can be derived from the topology of the corresponding 
Feynman graph G. Cutting I lines of a given connected /-loop graph such that it becomes a 
connected tree graph T defines a chord C(T,G) as being the set of lines not belonging to this 
tree. The Feynman parameters associated with each chord define a monomial of degree /. The 
set of all such trees (or 1 -trees) is denoted by T\. The 1 -trees T ET\ define U as being the sum 
over all monomials corresponding to the chords C (T, G) . Cutting one more line of a 1-tree leads 
to two disconnected trees (T\, Tz), or a 2-tree. T2 is the set of all such pairs. The corresponding 
chords define monomials of degree / + 1 . Each 2-tree of a graph corresponds to a cut defined 
by cutting the lines which connected the two now disconnected trees in the original graph. The 
square of the sum of momenta through the cut lines of one of the two disconnected trees T\ or T2 
defines a Lorentz invariant 

st = ( £ Pj) ■ (25) 

\jec(T,G) J 

The function Jo is the sum over all such monomials times minus the corresponding invariant. 
The function J is then given by Jo plus an additional piece involving the internal masses mj. In 
summary, the functions U and J are obtained from the graph as follows: 



« = e n 



X J 



j£C(T,G) 



'Jo 



(T u T 2 )e* 2 L ;ec(r!,G) 



1/ ' 



7 = To + uY^Xjm) . (26) 

7=1 
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In general, U is a positive semi-definite function. Its vanishing is related to the UV sub- 
divergences of the graph. Overall UV divergences, if present, will always be contained in the 
prefactor r(v — ID/ 2). In the Euclidean region, J is also a positive semi-definite function of 
the Feynman parameters Xj. The Euclidean region is defined as the region, where all invariants 
st are negative or zero. The vanishing of f is related to infrared divergences. Note that this is 
only a necessary but not sufficient condition for the occurrence of an infrared singularity. If an 
infrared singularity occurs or not depends in addition on the external kinematics. 

As an example for the functions U and J we consider the graph in fig. [U For simplicity we 
assume that all internal propagators are massless. Then the functions U and J read: 



Here we used the notation that Xy... r = xt +Xj + ... +x r . 

Before we close this section, let us consider one special case, where the following three 
conditions are met: 

1. The graph has no external lines or all invariants sj are zero. 

2. All internal masses mj are equal to /j. 

3. All propagators occur with power 1, i.e. Vj = 1 for all j. 
In this case the Feynman parameter integral reduces to 



This integral is a Igusa local zeta function (when viewed as a function of D/2) and has been 
studied by Belkale and Brosnan in [35]. 

4 Reduction to scalar integrals 

In the previous section we discussed scalar loop integrals, where the numerator of the integrand is 
independent of the loop momentum. This is not the most general case. A priori more complicated 
cases, where the loop momentum appears in the numerator might occur. However, there is a 
general reduction algorithm, which reduces these tensor integrals to scalar integrals [50,51]. The 
price we have to pay is that these scalar integrals involve higher powers of the propagators and/or 
have shifted dimensions. Therefore we considered in eq. (flTI) shifted dimensions and in eq. (fT8~l) 
arbitrary powers of the propagators. As a consequence, integrals of the form as in eq. ([TBI) are 



U = X\ 5 X23 + Xl5*46 +*23*46, 

7 = (xiX3X4+X 5 X 2 X 6 +XiX 5 X 2 346) (~Pl) 

+ (x 6 X3X 5 +X4X1X2 +X/PC6-X1235) {~Pl) 
+ (X2M*5 +X3XlX 6 +X2X3X U 56) (-pf) ■ 



(27) 




(28) 
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the most general loop integrals we have to consider. In this section we review the basic features 
of the reduction algorithm. 

We recall that the D-dimensional loop integrals over kj can be done loop by loop by Feynman 
parameterisation, completing the square in the denominator and a shift in the loop integration 
variables. The presence of the loop momentum in the numerator has the following effect: After 
the shift of the loop momentum, the (shifted) loop momentum as well as the Feynman parameter 
appear in the numerator. Let us first discuss the presence of the shifted (loop) momentum. We 
have to consider integrals of the form 

d D k 

*f f(k 2 ), (29) 



i% D l 2 



where f(k 2 ) is a function depending only on k 2 , but not on the individual components kF. Inte- 
grals with an odd number of the loop momentum in the numerator vanish by symmetry, while 
integrals with an even number of the loop momentum can be related by Lorentz invariance to the 
following integrals: 



dDk kWf(k 2 ) = -V [ ^(-k 2 )f(k 2 ), (30) 



/ 



/ iit D / 2 " " ' y " ' B° J in"/ 2 



The generalisation to arbitrary higher tensor structures is obvious. Recalling eq. (|2lJ) we observe 
that the dependency of the result on additional factors (— k 2 ) a occurs only in the combination 
D/2 + a, apart from a trivial factor T(D/2 + a) /T(D/2). Therefore adding a power of (—k 2 ) to 
the numerator is equivalent to consider the integral without this power in dimensions D + 2. 

In addition, shifting the loop momentum like in k' = k — xp introduces for tensor integrals 
the Feynman parameters Xj in the numerator. From the formula ([191 we observe that a Feynman 
parameter x in the numerator is equivalent to raising the power of the original propagator by one 
unit: V — ► v + 1 . Therefore we can relate an integral, where a Feynman parameter occurs in the 
numerator to a scalar integral, where the corresponding propagator is raised to a higher power. 

In summary, we can express all tensor integrals in terms of scalar integrals. The price we 
have to pay is that these scalar integrals involve higher powers of the propagators and/or have 
shifted dimensions. Therefore integrals of the form as in eq. (|22l) are the most general integrals 
which we have to consider. 



5 The main theorem 

Let us consider a general scalar multi-loop integral as in eq. (1221) 



Ig 



(a, 2 ) j d»*8(i-2>) i nv i JV _ lD/2 ■ (3D 

Xj>0 1 ' 
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Let m be an integer and set D = 2m — 2e. Then this integral has a Laurent series expansion in £ 



/ G = £ cj£j. (32) 

J=-v 

Theorem 1: In the case where 

1. all kinematical invariants st are negative or zero, 

2. all masses m, and /j are positive or zero (p ^ 0), 

3. all ratios of invariants and masses are rational, 

the coefficients Cj of the Laurent expansion are periods. 

Conditions (1) and (2) define the Euclidean region, where the integral does not possess an 
imaginary part. Condition (3) restricts all ratios to rational numbers. Since the rational numbers 
are dense in the real numbers, this is no severe restriction. 

Before we prove this theorem, we comment on two special cases. We first discuss the coef- 
ficient co in the case where the integral is convergent (and all propagators occur with power 1). 
In this case it is obvious, that cq is a period. An example is given by the one-loop three-point 
function with three external masses and vanishing internal masses. This integral is finite and no 
regularisation is needed. Therefore the Laurent series starts at e . The coefficient cq is given by 

f 3 8(1 -X\ -X2-X3) 



Xj >o 



d\X2X3 + Q2X-iX\ + CI3X1X2 



-(P2+P3) 2 ~{P3+Pl) 2 -(PI+P2) 2 

a\ = = , ci2 = ~ , «3 = = . (JJ) 

H L p z fj z 

By assumption, the parameters a\, 0,2 and are positive rational numbers. The integrand 
1/(01x2X3 +CI2X3X1 +CI3X1X2) is clearly a rational function, integrated over a semi- algebraic set. 
This shows that co is a period. 

The second case is given by integrals, where the conditions discussed at the end of section[3] 
are met: The graph has no external lines or all invariants sj are zero, all internal masses mj are 
equal to /j and all propagators occur with power 1 . In this case the integral reduces to a Igusa 
local zeta function of the form 

I G = f d n xh{\-j^Xi)U- D l 2 (34) 

Xj >0 i=1 

and it has been shown by Belkale and Brosnan [35] that the coefficients of the Laurent expansion 
are periods. 

We will actually prove a stronger version of theorem 1 . Consider the following integral 

j = r ^xsa-f^fn^^niPyW]^. (35) 
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The integration is over the standard simplex. The a's, b's, d's and f's are integers. The P's are 
polynomials in the variables x\, x n with rational coefficients. The polynomials are required 
to be non-zero inside the integration region, but may vanish on the boundaries of the integration 
region. To fix the sign, let us agree that all polynomials are positive inside the integration region. 
The integral J has a Laurent expansion 



Theorem 2: The coefficients Cj of the Laurent expansion of the integral J are periods. 

Theorem 1 follows then from theorem 2 as the special case a,- = V; — 1, bi = 0, r = 2, Pi = 11, 
p 2 = j , d\ + £/i = v - (/ + l)D/2 and d 2 + zfi = ID/2 - v. 

Proof of theorem 2: To prove the theorem we will give an algorithm which expresses each co- 
efficient Cj as a sum of absolutely convergent integrals over the unit hypercube with integrands, 
which are linear combinations of products of rational functions with logarithms of rational func- 
tions, all of them with rational coefficients. Let us denote this set of functions to which the 
integrands belong by M : Linear combinations of products of rational functions with logarithms 
of rational functions, all of them with rational coefficients. 

The unit hypercube is clearly a semi-algebraic set. From section |2] we know, that absolutely 
convergent integrals over semi-algebraic sets with integrands from the set M are periods. In ad- 
dition, the sum of periods is again a period. Therefore it is sufficient to express each coefficient 
Cj as a finite sum of absolutely convergent integrals over the unit hypercube with integrands from 
M . To do so, we use iterated sector decomposition [36-40]. The algorithm is described in detail 
in [36]. We proceed through the following steps: 

Step 0: Starting from the original integral J in eq. (|35T) we first convert all polynomials to homo- 
geneous polynomials. Due to the presence of the delta-function we have 



and we can multiply each term in each polynomial Pj by an appropriate power of x\ +x 2 + ■ ■ ■ +x n . 

Step 1: We can now assume that all polynomials are homogeneous. We then decompose the 
integral into n primary sectors as in 



J 



(36) 



j=jo 



1 



X\ +X 2 + ...+X n . 



(37) 




(38) 



In the Z-th primary sector we make the substitution 



xiXj for j 7^ 



(39) 
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and integrate out the variable x\ with the help of the delta-function. Each primary sector is 
now a (n — 1) -dimensional integral over the unit hypercube. Note that in the general case this 
decomposition introduces an additional polynomial factor 



n \ n r 

1 + L x j ' c = ~ n - E ( fl «' + zbi ) - E ( c j + ed j) » ( 4 °) 

where hj is the degree of the homogeneous polynomial Pj. After this step and after a relabelling 
we deal with integrals of the form 

l 

[d n xflxf +eb >fl[Pj(x)] d ^, (41) 

o !=1 j =i 

where the integral is now over the unit hypercube and the polynomials are positive semi-definite 
functions on the unit hypercube. Zeros may only occur on coordinate subspaces. Note that in 
general the polynomials Pj are no longer homogeneous. 



Step 2: We decompose the primary sectors iteratively into sub-sectors until each of the poly- 
nomials is of the form 

P = Cx^...x" n l "(l+P'(x)), (42) 

where P'(x) is a polynomial in the variables Xj not containing a constant term and C is a rational 
number. The term 1 + P'(x) is now a positive definite function on the unit hypercube. If P is 
of the form (|42|) . we say that P is monomialised. In this case the monomial prefactor x'" 1 ...x 1 "" 
can be factored out and the remainder contains a constant term. To convert P into the form (1421) 
we choose a subset S = {CL\, oc^} Q {1, ...n} according to a strategy. We decompose the 
fc-dimensional hypercube into k sub-sectors according to 

fd n x = t, [d n xfl Q(x a ,>x ai ). (43) 

In the Z-th sub-sector we make for each element of S the substitution 

Xai = x ai x a . fori^l. (44) 

This procedure is iterated, until all polynomials are of the form (l42l) . It is important to show 
that this can always be achieved in a finite number of iterations. In ref. [36] we gave a strategy 
for choosing the subset S, such that all polynomials are monomialised in a finite number of 
steps. This was done by relating the problem of monomialising the polynomials to the problem 
of the resolution of singularities of an algebraic variety over a field of characteristic zero. The 
monomialisation of a polynomial is equivalent to the special case, where the algebraic variety is 
defined through a single polynomial. For this special case Hironaka [41] invented his polyhedra 
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game to illustrate the challenge to find a constructive proof. Spivakovsky [42] was the first to 
give a winning strategy for Hironaka's polyhedra game and this winning strategy can be used 
in our case to ensure that the polynomials are monomialised in a finite number of steps. We 
review the relationship between Hironaka's polyhedra game and the problem of monomialising 
a polynomial in appendix lAl 

At the end of step 2 we have a finite number of sub- sector integrals. Each sub- sector integral 
is of the form as in eq. (|4TT) . where every Pj is now different from zero in the whole integration 
domain (including the boundaries). Hence the singular behaviour of the integral depends entirely 
on the factor 

f\x? +£bi . (45) 

Step 3: For every Xj with aj < we perform a Taylor expansion around xj = in order to extract 
the possible £-poles. If we consider for the moment only one parameter Xj we can write the 
corresponding integral as 



I 

where we defined l^ p > = d/dxjl(xj) 



(46) 



. The remainder term 



xj=0 

i( R \ Xj ) = /(*/)-£ (47) 

p=Q P- 

is by construction integrable in Xj and does not lead to £-poles from the xj -integration. The 
integration in the pole part can be carried out analytically: 

Idxx^il^ - ! — (48) 

J J J pi ~ aj + bje + p+1 p\ ■ (48) 
o 

This procedure is repeated for all variables xj for which aj < 0. At the end of step 3 we obtain a 
finite sum of integrals of the form 

l 

K(E) = -J- fd n xF(x,£), (49) 

g(e)J 

with 



F(x,t) = E/y(x,e), fj(x,z)=gj(E)Uxi II 
j=i i=i k=i 



(50) 
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Here, g(e) and gj(z) are polynomials in £ with integer coefficients. g(e) contains the explicit 
poles in e from the integration in eq. (|48T) . P[(x) is a polynomial with rational coefficients, 
non- vanishing on the unit hypercube. Further we have aj,bi,dj,fi E 7h. The integral in eq. (|49l) 

l 

f d n xF(x,E) (51) 
o 

is convergent by construction for all £ in a neighbourhood of £ = 0. In one variable this integral 
is of the form 

l 

J dxx £b R(x,e), (52) 
o 

where the function R(x,E) does not contain any singularities on the integration domain and is 
therefore bounded. Therefore the integral in eq. (|52|) is absolutely convergent for all £ with 

|e|<|l/fc|. 

Step 4: It remains to expand K(e), 1 /g(e) and F(x, e) in £: 

l 

oo -i oo to r p 

K(e)=Y,K r e r , -7-r = J>e r , F(x,e) = £F r (x)e r , K r = £ g s d n x F r _,(x) . (53) 

r=A #W r=A r=0 i=A ^ 

The expansion of the functions l/g(e) and gy(e) yields rational numbers, for the other terms we 
have 

k=0 

oo rk 

[P(x)] d+£ f = [P(x)} d Z J -(MP(x))) k e k . (54) 

k=o 

Therefore the integrand F r - S (x) belongs to the set M . 

The integrals over F r (x) are absolutely convergent: In each variable we have integrals of the 
form 

l 

j dx (lnx) k R r (x), fee N , (55) 
o 

where the function R r (x) does not contain any singularities on the integration domain and is 
therefore bounded. Therefore eq. (|55l) defines an absolutely convergent integral. This completes 
the proof. 
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6 Conclusions 



In this article we considered the Laurent expansion of multi-loop Feynman integrals in dimen- 
sional regularisation. We studied the integral in the Euclidean region and where all ratios of 
invariants and masses have rational values. We showed that in this case all coefficients of the 
Laurent series are periods (theorem 1). This is an important result for the theory of Feynman 
integrals, as it restricts the class of functions which can appear in the computation of Feynman 
integrals. The proof of theorem 1 follows from a stronger result, stated in theorem 2: Integrals 
of the form as in eq. (1351) have a Laurent expansion, where all coefficients are periods. The proof 
of theorem 2 is constructive and based on iterated sector decomposition. 

A Hironaka's polyhedra game 

In this appendix we show the equivalence of Hironaka's polyhedra game with the problem of 
monomialising a polynomial. 

Hironaka's polyhedra game is played by two players, A and B. They are given a finite set M 
of points m = (mi, m n ) in N" , the first quadrant of N B . We denote by A C W]_ the positive 
convex hull of the set M. It is given by the convex hull of the set 

(J (m + R n + ). (56) 

The two players compete in the following game: 

1. Player A chooses a non-empty subset S C {1, n}. 

2. Player B chooses one element i out of this subset S. 

Then, according to the choices of the players, the components of all (mi , . . . , m n ) G M are replaced 
by new points (m\, m' n ), given by: 

m'j = ntj, if j ^ z, 

where for the moment we set c = 1. This defines the set M'. One then sets M = M' and goes 
back to step 1 . Player A wins the game if, after a finite number of moves, the polyhedron A is of 
the form 

A = m + R n + , (58) 

i.e. generated by one point. If this never occurs, player B has won. The challenge of the polyhe- 
dra game is to show that player A always has a winning strategy, no matter how player B chooses 
his moves. 
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Let us discuss the relation of Hironaka's polyhedra game to the sector decomposition of 
multi-loop integrals. Without loss of generality we can assume that we have just one polynomial 
P in eq. (1351) . (If there are several polynomials, we obtain a single polynomial by multiply- 
ing them together. As only the zero-sets of the polynomials are relevant, the exponents can be 
neglected.) The polynomial P has the form 

f=l 

The n-tuple m$ = (nif\ m$ j defines a point in N" andM = jm^ 1 ), ...m^| is the set of all 

such points. Substituting the parameters xj according to equation (|44|) and factoring out a term xf 
yields the same polynomial as replacing the powers mj according to equation (1571) . In this sense, 
one iteration of the sector decomposition corresponds to one move in Hironaka's game. Reducing 
P to the form ((42]) is equivalent to achieving (1581) in the polyhedra game. Finding a strategy which 
guarantees termination of the iterated sector decomposition corresponds to a winning strategy for 
player A in the polyhedra game. Note that we really need a strategy that guarantees player A's 
victory for every choice player B can take, because the sector decomposition has to be carried 
out in every appearing sector. In other words, we sample over all possible decisions of B. 

B The prefactors 

In this appendix we discuss the prefactors, which we dropped in eq. (PT8T) . From the Feynman 
rules we obtain for a loop integral 

with v = Vi + ... + V n . Here we included for each loop an additional coupling factor g 2 relative 
to a corresponding tree graph and a factor (e^ / 4 / n) £ related to the MS-scheme [52]. We have 

h = C G I G: (61) 

where Iq is defined by eq. (TT8T) and the prefactor Cq is given by 

Cc = ^yr^r^^Y^zim. (62) 

y4Kj nr(v ; ) 



With D = 2m — 2e we obtain 



Cg = T^S^ Er(v "' m+ ' e)) ^' (63) 
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Since all Vj £ N, the term (-l) v i n+ '/nr(Vj) is obviously a period. We can further show that 
the second term has a Laurent expansion in £, where all coefficients are periods: 

^ e r(v-/m+/ e ) = ^r(i+te)- 5^±^ . (64) 

Using the functional equation r(x-|- 1) = xT(x) one shows that T(v — lm + /e)/T(l + le) is a 
rational function in £ with rational coefficients. The expansion of T(l + e) is given by 



n=2 H 




r(l+e) = exp(-y £ £+£^-^^£"), (65) 

and therefore 

e z ^ e r(l+/£) = exp ( V ^-^ n e n ) . (66) 




The zeta values C, n are periods, therefore the product e ly£Z T(v — lm + le) has a Laurent expansion 
in £, where all coefficients are periods. Note that Euler's constant Je, which is conjectured not to 
be a period, drops out in this expression. 

This leaves the factor g 21 / (4n) lm . It is not known whether l/% is a period or not. Kontsevich 
and Zagier [44] consider therefore the extended period ring P = P[1/tt], obtained by adding 
1/71 to the ring of periods. In any case it is common practice to quote results of a perturbative 
calculation as an expansion in g 2 / (An)" 1 : For D = 4 — 2e and with a = g 2 / (4%) the result is often 
expressed as 



o 



a „ ( a \2 



(67) 



In the coefficients Cj the factors g 2 / (4tc) 2 = a/ (47l) have been explicitly factored out. The results 
of this paper apply to the coefficients Cj. 



C Other algorithms for the sequence of blow-ups 

Triggered by the response of a referee we show in this appendix that the algorithm for the con- 
struction of the sequence of blow-ups given by Bloch, Esnault and Kreimer in [43] does in gen- 
eral not monomialize a polynomial. We do this by giving an explicit counter-example. Consider 
the polynomial 

p = X2+X1X3 (68) 

with homogeneous coordinates x\, X2, X3, X4 on P 3 . We will examine the polynomial P in a 
neighbourhood of (xi,;t2>*3>*4) = (0,0,0, 1). Let us denote coordinate sub-spaces by L ;J : X[ = 
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Xj = and : xt = Xj = Xk = 0. The algorithm of [43] determines first the set S of coordinate 
subspaces, on which P vanishes and which are maximal. In our case 

S = {Li2,L 23 }. (69) 

One then forms all possible intersections of elements of S, which defines the set 7 : 

7 = {L 12 ,L 23 ,L 123 }. (70) 

7min contains the minimal elements of 7 '■ 

7,nin = {L m }- (VI) 

The elements of 7min define the centres of the first round of blow-ups. In our case there is only 
one centre Li 23 . The strict transforms of the elements of 7 which are not in 7 m i n define 7\. By 
abuse of notation we denote the strict transforms of Ly also by L, ; -: 

7\ = 7\7min = {Ln,L 2 3}. (72) 

The minimal elements of 7\ define the centres of the second round of blow-ups. In our case 
7\min — 7\- This procedure is continued by defining 7i = 7\\7\,min etc - untl l 7j — ®- I n 
our case 72 — and we stop after the second round of blow-ups. Let us now investigate how 
the polynomial transforms under the blow-ups: The first one blows up the centre Li 23 . In a 
coordinate patch where x\ > * 2 and x\ > jc 3 we have with x\ 

p = x j + xix 3 = x\ 2 (x' 2 2 + x' 3 ) = x[ 2 P' . (73) 

We then have to blow up the centres of 7\min- The strict transform of L\2 is not contained in the 
above coordinate patch, but the strict transform of L 23 is. Let us now consider what happens if 
we blow up the centre given by the strict transform of L 23 . Looking in the coordinate patch with 
x'2 > x'3 we have with x[ = x'(, x' 2 = x' 2 \ x' 3 = x'^x": 

p' = 4(4+4) =4p". (74) 

We see that P" is not in monomialized form. The proof of proposition 7.3 in [43] relies crucially 
on the assumptions (i) that only integrals of the form as in eq. (1281) are considered, (ii) that the 
graph under consideration is a graph with n loops and In edges, (iii) that all proper subgraphs 
are convergent. In the present paper we do not make these assumptions. In fact in our theorem 2 
arbitrary polynomials may appear, which need not be related to Feynman graphs. Our theorem is 
therefore a much stronger statement. In particular polynomials like eq. (1681) are allowed, which 
vanish along one direction quadratically and along other directions linearly. 
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